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Abstract — In this work, we present an analytical method 
for calculating the axial and radial components of the 
magnetic field created by a permanent magnet in the form 
of ring used for magnetron Radio Frequency (RF) reactor, 
the calculation is based on the complete elliptic integral of 
first kind and also heuman lambda function by using 
Linux 64-bit. This method allows us to accurately assess 
the magnetic field, and also allows us to reduce the 
computational time. 
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I. INTRODUCTION 


The magnetrons Radio Frequency (RF) discharges 
are used in particular for the deposition of thin metal or 
insulating layers for nano-components used in the field 
of micro-electronics, for surface treatment, and also, for 
the Deep Reactive Ion Etching (DRIE Technology) of 
the Silicon to manufacture capacitors with high density 
and the depositing thick insulation or thick metal to 
manufacture micro-coils. The references [1-6] are 
examples of theoretical and experimental work in this 
area. In these reactors, a permanent magnet is placed 
next to the cathode surface in order to make the 
electrons stay longer in the plasma, and also to confine 
these electrons before they are collected or recombined 
at the electrodes. This phenomenon allows increasing 
the plasma density [2] which is much lower in the 
absence of magnetic field due to the very low gas 
densities [3]. Under these conditions, the ions generated 
in the plasma effectively contribute to the cathodic 
sputtering of target and forming in the gas phase the 
necessary precursors to deposit the thin films on the 
substrate placed on the anode electrode. The 
optimization of such a deposit can be made using 
experimental diagnostic tools supplemented with tools 
for modeling the dynamics of the discharge and the 
kinetics of the non-equilibrium plasma formed. These 
models require absolutely a priori knowledge of the 
magnetic field that exists within the plasma reactor. 
That is why the purpose of this paper is to propose an 
analytical model for calculating the magnetic field 
created by one a permanent magnet ring radially 
magnetized to be used in the future in a particle model 
for the simulation of a magnetron reactor. Several 
analytical and numerical approaches have been made to 
the calculation of the magnetic field created by the 
magnets [7-11]. 


The main objective of this work is to calculate the 
axial and radial components of the magnetic field in the 
volume of the discharge created by a permanent magnet 
in a ring shape radially magnetized. This calculation is 
based on the complete elliptic integral of first kind also 
heuman lambda function. Figure | shows a scheme of a 
magnetron RF reactor, the domain of the calculation is 
defined by the volume of the reactor constituted by two 
parallel and circular electrodes of radius R, and the 
distance between them is H,, the distance between the 
cathode and the permanent magnet is (H-h) where h is 
the amplitude of the magnet (see figure.2). 
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Fig. 1 Scheme of a magnetron RF reactor 
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Fig.2. Scheme of design of the magnet with the reactor 


II. BASIC EXPRESSIONS OF COMPUTING THE MAGNETIC 
FIELD 


The calculation of the magnetic field created by a 
permanent magnet can be obtained by using the 
Coulomb method current model [8], [11]. In this paper 
we used the coulomb method, in the case of a radial 
magnetization (see fig.2), the magnetic field at any 
point of coordinates M (r,z) created by the inner face of 
magnetization +o and outer face of magnetization -o 
respectively is given by the following expressions: 


= 2 yh 6-22 Wa l 
B*(r,z) = = ee ey ) 
ae SS 3 1 1 
4a z=0 0=0 | 
: = z=h O=20 ava 2 
Bo (r,z) = = i) i) ae r,.dz,.d0 @) 
4n z=0 6=0 |P a | 
B (r,z) = B* (r,z)+ Bo (r,z) (3) 


PM = (r-r,cos@)i,—1r,sin0i7,+(z-2,)i, (4) 


Where P(7r),Z1,8 ) is point located on the inner or 
outer face of the permanent magnet defined by its radial 
position 7; (r=Ri, for the inner face and r,;=R,, for the 
outer face), axial position z, and polar angle 6 and the 
equation (1) can be re-written as follows: 
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In this equation, the first term represents the radial 
component of the magnetic field created by the inner 
face of the permanent magnet B,, the second term 
represents the tangential component Bo’ and third term 
represents the axial component B,. Due to the 
symmetry, the tangential component is zero (By = 0). 
The same steps may be followed in order to find the 
radial component (B8,) and the axial component (B,) 
created by the magnetization of the outer face of the 
magnet ring. By substituting the variable 0 = 229 in 
the equation (5) and after some mathematical 
calculation we obtain the following expressions: 
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K(k) is the complete elliptic integral of first kind and 
Ao(f,k) is the heuman lambda function where: 


K (k)=R, (0,1-?,1) (8) 
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Ry and R; are respectively the symmetric integrals 
of first kind and third kind where: 
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Where the square root is taken real and positive if x, 
y, Z are positive and varies continuously when x, y, z 
become complex. The integral is well defined if x, y, z 
lie in the complex plane cut along the non positive real 
axis (henceforth called the “cut plane’), with the 
exception that at most one of x, y, z may be 0. The same 
requirements apply to the symmetric integral of the 
third kind, where p # 0 and the Cauchy principal value 
is to be taken if p is real and negative. A degenerate 
case of R; that embraces the inverse circular and 
inverse hyperbolic functions. 

The reader can find the algorithm of these kinds of 
integrals in the reference (12). For accurate results, 
especially at the vicinity of singular points, the code of 
calculation was performed under linux 64-bit. 


III. SIMULATION RESULTS 


The dimensions of the geometry used in this 
calculation are as follows, H = 3.5 cm, h = 3 cm, R, = 6 
em, H, = 2.5 cm, Rj, = 2 cm and R,, = 4 cm. The 
magnetization is considered constant o = IT for the 
inner face and o = -IT for the outer face (IT = 10° 
Gauss). 

In figure 3 and figure 4 respectively, we show the 
contours of radial component and axial component of 
magnetic field created by the permanent magnet ring in 
all volume shown in figure 2 and figure 5 shows the 
contours of the total magnetic field in the same volume. 
It is clear that the points which are located in the 
corners of the permanent magnet ring are singular 
points. These points are treated by using the Hospital 
rule. Figure 6 illustrate the spatial variation of the 
radial component of the magnetic field in the volume of 
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reactor and created by the permanent magnet ring 
radially magnetized, this component is calculated only 
in the bulk of reactor which located by (H-h) = 0.5 cm 
from the magnet ring in the upper side. This component 
is large in the vicinity of the magnetic ring and the 
location of the target. It is necessary to say that this 
component has the major effect on plasma confinement 
by making the particles charged stay longer in this 
region and as result of that, it is effectively contribute to 
the cathodic sputtering. Figure 7 shows the spatial 
variation of the axial component of the magnetic field 
in the volume of reactor and created by the permanent 
magnet ring radially magnetized, this component does 
not have an effect on the dynamic of plasma. 


Finally, figure 8 shows the spatial variation of the 0 1 2 2 4 5 6 
total magnetic field in the bulk of the asymmetric r (cm) 
magnetron reactor, it represents the magnitude of the 
total magnetic field. Fig.5. Contours of the total magnetic field created by 


permanent magnet ring radially magnetized 
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Fig.3. Radial component contours of magnetic field created Fig.6. Spatial variation of the radial component of the 
by permanent magnet ring radially magnetized magnetic field in the bulk of reactor 
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Fig.4. Axial component contours of magnetic field created Fig.7. Spatial variation of the axial component of the 
by permanent magnet ring radially magnetized magnetic field in the bulk of reactor 
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Fig.8. Spatial variation of the total magnetic field in the 
bulk of reactor 


IV. CONCLUSION 


In this paper we presented the analytical expressions 
for the calculation of the components of the magnetic 
field created by one permanent magnet ring, this 
magnet can be used in the magnetron reactor for the 
confinement of the plasma. These expressions are valid 
to calculate the magnetic field at any regular point in 
the space. This kind of calculation is simple but it is 
faster in comparison with the usual numerical 
calculations of magnetic field. 
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